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Mathematical  Modeling  and  Simulation  of  Flexible  Structures 


Abstract 

The  present  documment  summarizes  the  main  results  obtained  during  the  period 
1988/1990  in  areas  related  to  the  mathematical  modeling  and  numerical  simulation  of  the 
dynamics  of  flexible  structures  undergoing  large  overall  motions.  This  is  the  final  report 
of  a  two  year  research  effort  sponsored  by  AFOSR  under  current  grant  DJA/AFOSR  89- 
0294  with  Stanford  University,  as  a  continuation  of  the  effort  initiated  in  1986  under  greint 
DJA/AFOSR  86-0292  with  Stanford  University. 

The  work  reported  on  herein  is  part  of  a  multidisciplinary  effort  aimed  at  the  de¬ 
velopment  of  analytical  and  computational  tools  for  the  modeling  ?md  simulation  of  very 
flexible  structures.  This  research  is  concerned  with  three  main  areas:  (i)  Mathematical 
sispects.  with  special  e-  phasis  on  nonlinear  stability  analysis;  (ii)  Modeling  aspects,  with 
emphasis  in  very  flex.^ie  slender  structures;  and  (iii)  Numerical  analysis  aspects,  with 
emphasis  on  methodologies  suitable  for  large  scale  and  accurate  simulation.  The  present 
research  effort  has  resulted  in  36  publications  (see  Appendix  I),  two  conferences,  and  a 
strong  collaboration  (in  a  related  research  effort  also  sponsored  by  AFOSR)  with  groups 
UC  Berkeley  and  U.  Maryland. 
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1.  Introduction 

The  present  report  gives  an  overview  of  the  main  results  obtained  during  the  past  two 
years  of  research  sponsored  by  AFOSR  imder  current  grant  DJA/AFOSR  89-0294  with 
Stanford  University.  This  effort  is  a  continuation  of  a  previous  two-yeax  research  effort 
also  sponsored  by  AFOSR  imder  grant  DJA/AFOSR  86-0292  with  Stanford  University. 
The  goal  of  this  work  is  the  development  of  analytical  and  computational  tools  for  design, 
analysis  eind  simulation  of  the  dynamics  of  very  flexible  structiures  imdergoing  large  overall 
motions. 

The  organization  of  this  report  is  as  follows.  Section  2  gives  a  brief  outline  of  the 
scope  of  the  proposed  research  along  with  a  summary  of  the  research  goals  set  forth  in  our 
proposal  to  AFOSR.  Section  3  gives  a  summary  of  the  most  important  results  obtained  in 
the  course  of  this  investigation  which  fall  within  three  main  areas;  (i)  Mathematical  aspects 
involved  in  the  stability  analysis  of  the  mechanical  systems  of  interest,  (ii)  Mathematiced 
modeling  of  the  physical  systems  of  interest;  in  particular,  very  flexible  structures,  and  (iii) 
Numerical  analysis  and  computational  aspects  involved  in  the  simulation  of  the  dynamic 
response  of  these  systems.  Finally,  Section  3  gives  a  brief  outline  of  some  lines  of  research 
to  be  pm-sued  in  a  continuation  of  this  effort.  Further  details  on  the  scope  of  this  proposed 
research  are  given  in  a  follow-up  proposal.  Appendix  I  contains  a  complete  relation  of 
the  36  publication  produced  under  the  sponsorship  of  AFOSR  since  the  initiation  of  this 
research. 


2.  Scope  of  the  Proposed  Research 

As  stated  in  om  proposal  to  AFOSR,  the  centred  theme  of  the  present  investigation 
is  the  formulation,  mathematical  and  numerical  einalysis  of  models  for  flexible  multi-body 
systems,  suitable  for  precision  control  and  amenable  to  large-scale  simidation.  In  particu- 
leur,  we  are  concerned  with  so-called  geometrically  exact  structmal  models,  a  terminology 
introduced  by  s.  ANTMAN  and  J.  SIMO,  for  which  no  restriction  is  placed  on  the  degree 
of  allowable  flexibility  in  the  structure.  From  an  engineering  point  of  view,  this  class  of 
models  allow  the  correct  modeling  of  important  phenomena  such  as  large  deformations, 
nonlinear  vibration,  bifmcations  and  structural  daunage  which  may  arise,  for  instance,  as 
a  result  of  control  system  breakdown.  We  remark  that  geometrically  nonlinear  effects  play 
an  important  role  in  a  number  of  engineering  applications,  particularly  those  involving  ro¬ 
tating  structmes,  and  their  prediction  requires  properly  formulated  nonlinem  models.  In 
fact,  as  demonstrated  in  SIMO  &  VU-<JUOC  [1988b],  use  of  linearised  theories  can  produce 
completely  misleading  and  erroneous  results. 

Within  the  preceding  context,  the  following  three  basic  areas  of  investigation  were 
proposed  for  the  period  1938/1990  (see  Section  5  of  our  renewzJ  proposal); 

i.  Formulation  and  analysis  of  geometrically  exact  shell  models.  In  particulEu-  two 
concrete  objectives  were  proposed  within  this  area  of  research.  First,  the  formulation  of 
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geometrically  exact  shell  models  which  account  for  thickness  change.  This  class  of  mod¬ 
els  axe  of  considerable  interest  in  technological  applications  concerned,  for  instance,  with 
composite  shell  structures  where  the  prediction  important  effects  such  as  delamination  re¬ 
quire  the  use  of  these  enheinced  models.  Second,  the  formulation,  algorithmic  development 
and  numerical  analysis  aspects  involved  in  the  implementation  of  inelaisticity  and  damage 
constitutive  models  for  geometrically  exact  shells. 

ii.  Numerical  analysis  of  the  dynamics  of  shell  models.  It  was  proposed  to  address  in 
detail  numerical  analysis  aspects  involved  in  the  actual  implementation  of  geometrically  ex¬ 
act  shell  models.  In  particular,  the  development  of  exact  and  singularity-free  configuration 
update  procedures,  the  investigation  of  appropriate  finite  element  discretizations  suitable 
for  effective  numerical  simulations,  and  the  use  of  covariant  linearization  techniques  to 
derive  exact  tangent  operators  which  plays  a  crucial  role  in  incremental,  iterative,  solu¬ 
tion  strategies.  The  overall  goal  here  is  the  analysis  and  development  of  computational 
capabilities  suitable  for  large-scale  simulation. 

Hi.  Stability  Analysis.  The  Energy- Momentum  method.  The  two  areas  alluded  to 
above  are  concerned  primarily  with  the  development  of  mathematical  models,  numerical 
analysis  techniques  and  computational  tools  for  the  integrated  design  of  flexible  structures. 
In  addition,  it  was  proposed  to  continue  the  efforts  on  frmdamental  areas  of  hamiltonian 
and  continuum  mechanics,  with  emphasis  on  applications  to  stability  analysis  of  relative 
equilibria.  This  problem  is  of  particiilar  importance  in  areas  such  8is  spacecraft  dynamics 
where  effects  related  to  the  flexibility  of  the  appendages  (antenae,  solar  panels,  etc...,)  play 
a  critical  role  in  the  stability  of  the  system. 

Om  proposal  contained  a  fourth  area  of  research  concerned  with  the  formulation 
and  numerical  zinalysis  of  conserving  (symplectic)  time  stepping-algorithms  for  nonlinear 
mecheinics.  Recently,  we  have  obtained  new  and  rather  exciting  results  in  this  area,  which 
will  be  reported  on  in  the  follow-up  renewal  proposal. 


3.  Summary  of  research  objectives  accomplished 

During  the  last  year  we  have  made  considerable  progress  in  the  three  areas  of  research 
outlined  above.  The  main  developments  are  summarized  below  along  with  a  comment  on 
selected  publications  recorded  in  Appendix  I. 

i.  Formulation  and  analysis  of  geometrically  exact  shell  models.  The  two  objectives 
set  forth  imder  this  line  of  reseeu-ch  have  been  completed.  In  SIMO,  FOX  &:  RIFAI  [1990b], 
we  have  presented  and  developed  in  detail  a  fully  nonlinear  shell  theory  that  incorporates 
finite  thickness  changes.  A  comprehensive  modem  treatment  and  summary  of  the  current 
status  of  nonlinear  shell  theory  is  given  in  the  thesis  of  FOX  [1990].  In  SIMO,  MARSDEN 
&  KRISHNAPRASAD  [1988],  we  have  completed  our  program  on  the  systematic  study  of 
the  hamiltonian  structure  of  nonlinear  elasticity,  rods  and  shells.  These  results  play  a 
crucied  role  in  the  formulation  and  analysis  of  the  algorithms  described  below.  In  SIMO 
&  KENNEDY  [1990],  we  have  imdertaken  the  formulation,  numerical  amalysis  and  detailed 
implementation  of  inelastic  constitutive  models  formulated  in  stress  resultants.  This  paper 
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also  contains  a  nxunber  representative  large  scale  numerical  simulations;  for  instance,  the 
the  full  simulation  of  the  inelastic  collapse  of  cylindrical  shells  in  the  plastic  regime. 

ii.  Numerical  analysis  of  the  dynamics  of  shell  models.  The  objectives  and  program 
of  research  set  forth  vmder  this  line  of  investigation  have  also  been  completed.  In  paxticulax, 
in  SIMO,  FOX  h  RlFAl  [1990a, b],  we  have  provided  a  complete  numerical  ainalysis  and  algo¬ 
rithmic  treatment  of  nonlineax  (geometrically  exact)  shell  theory  incorporating  thickness 
changes.  This  includes  the  construction  of  singularity-free  shell  parametrizations  and  the 
development  of  effective  configuration  update  procedures  which  remain  exact  regardless 
the  magnitude  of  the  incremental  displacements  and/or  rotations.  Such  a  development 
exploits  in  a  crucial  manner  the  geometric  framework  developed  in  item  i  above,  along 
with  ideas  of  modern  differential  geometry.  Related  developments  for  nonlineax  rods  axe 
addressed  in  SIMO  &  VU-QUOC  [1988a,b,1990]. 

Simulations  of  the  dynamics  of  fiexible  structures  rely  crucially  on  finite  element 
discretizations  of  the  continuum  problem.  For  plates  and  shell  structures,  standard  finite 
element  methods  often  lead  to  discretizations  that  experience  severe  convergence  and  ac- 
cmacy  problems  in  the  thin  shell  limit.  In  the  past,  ad-hoc  methods  have  been  developed, 
typically  from  an  engineering  perspective,  to  circiimvent  these  well-known  difficvilties.  In 
SIMO  &  RIFAI  [1990]  a  general  class  of  stable  mixed  finite  element  methods  is  developed 
for  a  number  of  constrained  problems]  in  paxticulax,  very  thin  shell,  which  exhibit  excellent 
accuracy  characteristics  for  coaxse  meshes.  These  and  related  finite  element  techniques 
axe  employed  and  further  developed  in  the  specific  context  of  shells  in  SIMO,  FOX  &  RI¬ 
FAI  [1990a,b].  Further  applications  axe  considered  in  references  [28-32]  of  Appendix  I.  An 
overview  of  the  state  of  the  art  in  mathematical  and  numerical  analysis  aspects  related  to 
rods,  plates  and  shells  is  contained  in  the  recent  book  quoted  in  reference  [34]. 

iii.  Stability  analysis.  The  Energy-Momentum  method.  The  progress  made  in  this 
area  of  reseaxch  over  the  past  two  years  has  led  to  what  appears  to  be  significant  new  re- 
stolts  in  methodologies  for  the  stability  analysis  ( and  eventual  bifurcation)  of  hamiltonian 
systems  with  symmetry.  Initially,  it  was  proposed  to  consider  as  a  representative  model 
problem  the  stability  of  coupled  rigid  bodies/flexible  attachments  by  means  of  a  methodol¬ 
ogy  going  back  to  ARNOLD  [1966]  and  currently  referred  to  as  the  energy-Casimir  method. 
Remarkably,  however,  this  well-established  method  fails  for  realistic  models  of  the  flex¬ 
ible  appendage.  The  crucied  difficxilty  lies  precisely  in  the  actual  characterization  of  the 
Casimir  functions  associated  with  the  mathematical  model  of  the  appendage..  We  recall 
that  Ceisimirs  are  conserved  quantities  in  the  reduced  space;  e.g.,  for  rigid  bodies  the 
norm  of  the  total  angular  momentum.  For  nonlineax  elasticity,  rods  and  shells  no  explicit 
Casimirs  appear  to  be  known. 

In  SIMO,  POSBERGH  &  MARSDEN  [1990a],  the  difficulties  alluded  to  above  axe  com¬ 
pletely  bypassed  by  exploiting  directly  the  veuiational  chareicterization  of  relative  equilib¬ 
ria  as  critical  points  of  the  energy-momentum  map,  and  working  directly  on  the  canonical 
phase  space.  Casimir  functions,  which  as  pointed  out  above  may  difficult  to  chaxEicterize 
or  may  not  exist  at  all,  play  no  role  in  this  approach.  The  methodology  is  further  refined, 
generalized  and  successfully  applied  to  specific  examples  in  LEWIS  &  SIMO  [1990a].  Addi¬ 
tional  applications  and  mathematical  aspects  of  the  method  axe  considered  in  references 
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[22-26]  of  Appendix  I.  The  proceedings  of  the  1998  summer  conference  (see  reference  [35]) 
cont^un  a  number  of  related  papers  and  give  a  recent  overview  of  work  in  this  area. 

We  conclude  this  brief  overview  by  noting  the  positive  effect  on  the  research  reported 
on  above  derived  from  our  continued  interaction  with  J.E.  MARSDEN  at  Berkeley,  and  P.S. 
KRISHNAPRASAD  at  Maryland  in  a  related  effort  also  sponsored  by  AFOSR. 


3.  Proposed  areas  of  research  in  a  continuation  effort 

The  objective  of  this  Section  is  to  advance  a  number  of  lines  of  research  to  be  pursued 
in  a  continuation  of  this  investigation.  A  detailed  discussion  of  these  research  objectives 
is  deferred  to  a  follow-up  proposal  to  be  submitted  to  AFOSR  shortly  after  this  report. 
In  keeping  with  the  structure  of  this  report,  the  proposed  research  is  organized  along  the 
same  lines  considered  in  Section  2. 

i.  Mathematical  models  of  flexible  structure.  Realistic  models  of  spacecrafts,  space 
platforms  and  related  devices  invariably  involve  non-smooth  intersections  of  structural 
elements  such  as  plates,  rods  and  shells.  The  appropriate  mathematical  and  nvunerical 
modeling  of  these  intersection  problems  has  been  a  subject  of  considerable  interest  and 
some  controversy  in  recent  years.  In  paxtictdar,  shell  models  with  drill  rotations  have  been 
proposed  as  a  rather  practical  means  of  resolving  shell  intersection  problems.  We  plan  to 
investigate  the  mathematical  foxmdations  and  numerical  implications  of  this  class  of  widely 
used  structural  models. 

ii.  Numerical  analysis.  As  noted  in  our  proposal,  a  large  body  of  computational 
experience  indicates  that  standard  time-stepping  algorithms  are  not  suitable  for  the  long¬ 
term  simulation  on  nonlinear  dynamical  systems.  In  fact,  oiu:  recent  results  suggest  that 
non  of  the  existing  and  widely  used  time  integrators  in  structural  dynamics  conserve  fun¬ 
damental  constemts  of  motion  such  eis  angular  momentum  and  energy.  Recently,  we  have 
been  successful  in  developing  algorithms  for  rigid  body  dynamics  which  conserve  exactly 
these  constants  of  motion  without  resorting  to  Lagrange  multiplier  and  related  (awkward) 
devices.  Aside  from  its  practical  implications,  conservation  properties  lead  naturailly  to 
nonlinejir  notions  of  algorithmic  stability.  Our  objective  is  to  further  analyze  these  meth¬ 
ods  and  extend  them  to  more  sophisticated  and  mechemically  more  interesting  models  such 
as  rods  and  shells. 

Standard  iterative  solution  procedures  for  the  numerical  simulation  of  nonlinear  prob¬ 
lems  break  down  in  the  presence  of  nonlinear  effects  such  as  transfer  of  stability  and  bi¬ 
furcation  as  a  result  of  symmetry  breaking.  In  such  situations  which  arise,  in  particular, 
in  numerical  simulations  of  nonlinear  shells,  it  is  necessary  to  use  continuation  methods. 
Our  preliminary  work  in  WRIGGERS  &  SIMO  [1990]  suggests  that  it  is  possible  to  develop 
algorithms  of  this  type  which  allow  the  direct  computation  of  critical  points  by  iterative 
methods  that  exhibit  quadratic  rate  of  convergence. 

iii.  Geometric  methods  in  stability  and  bifurcation  analysis.  As  pointed  out  above, 
the  energy-momentum  method  developed  in  SIMO,  POSBERGH  &  MARSDEN  [1990]  offers 


6 


J.C.  Simo 


7 


a  general  methodology  for  the  systematic  stability  analysis  of  relative  equilibria  in  hamil- 
tonian  systems.  The  crucial  tool  in  this  approach  is  a  block-diagonalization  result  which 
decouples,  as  far  as  the  stability  analysis  is  concerned,  rotational  from  internal  or  vibra¬ 
tional  modes.  Note  that  these  modes  are  dynamically  coupled.  Nevertheless,  we  show  that 
a  change  of  variables  can  be  explicitly  constructed  such  that  the  second  variation  of  the 
energy-momentum  map  block-diagonalizes.  As  a  result,  the  stability  conditions  associated 
with  the  rigid  modes  become  explicit. 

The  preceding  block  diagonalization  result  has  important  implications  in  stability 
and  bifurcation  analyses  of  rotating  flexible  structures  that  we  plan  to  explore  in  detail 
in  the  future.  In  particular,  oxir  preliminary  work  in  references  [1]  and  [2]  suggests  a 
further  simplification  of  the  energy-momentum  method  crucial  in  concrete  applications. 
Furthermore,  the  coordinate  change  involved  in  the  block-diagonalization  theorem  has  the 
remarkable  property  of  providing  an  explicit  normal  form  for  the  linearized  dynamics  at  a 
critical  point.  Preliminary  results  indicate  that  this  property  plays  a  key  role  in  bifurcation 
analyses  at  a  relative  equilibrium. 


Appendix  I:  Publications  under  AFOSR  Support 


I.l.  Publications  in  Referee  Journals 

[1]  Simo,  J.C.,  D.  Lewis,  and  J.E.  Marsden,  [1990]  “Stability  of  relative  equilibria.  Part 
I;  The  reduced  energy-momentum  method,”  SUDAM  Report  No.  89-3.  Submitted 
to  Archive  for  Rational  Mechanics  and  Analysis. 

[2]  Simo,  J.C.,  T.  Posbergh,  and  J.E.  Marsden,  [1990a]  “Stability  of  relative  equilibria. 
Part  II:  Application  to  nonlinear  elasticity,”  SUDAM  Report  No.  89-4.  Submitted  to 
Archive  for  Rational  Mechanics  and  Analysis. 

[3]  Simo,  J.C.,  T.  Posbergh,  and  J.E.  Marsden,  [1990b],  “Stability  analysis  of  coupled 
rigid  bodies  and  geometrically  exact  rods;  Block  diagonalization  and  the  energy- 
momentum  method,”  Physics  Reports,  to  appear. 

[4]  Lewis,  D.,  and  J.C.  Simo,  [1990],  “Nonlinear  stability  of  pseudo-rigid  bodies,”  Pro¬ 
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